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29 ABSTRACT 

30 We present the light curves and spectral data of two exceptionally luminous gamma-ray 
outburts observed by the Large Area Telescope (LAT) experiment on board Fermi Gamma-ray 
Space Telescope from 3C 273 in September 2009. During these flares, having a duration of a 
few days, the source reached its highest 7-ray flux ever measured. This allowed us to study 
in some details their spectral and temporal structures. The rise and decay are asymmetric on 
timescales of 6 hours, and the spectral index was significantly harder during the flares than 
during the preceding 1 1 months. We also found that short, very intense flares put out the same 
time-integrated energy as long, less intense flares like that observed in August 2009. 

31 Subject headings: Galaxies: active — quasars: Individual: 3C 273 — Gamma rays: observa- 

32 tions — 



33 1. Introduction 

34 The LAT experiment, on board the Fermi Gamma-ray Space Telescope satellite, observes the entire 

35 sky, in the 0.02 - >300 GeV band, once every ~3 hours. It is providing the first collection of well sampled 

36 gamma-ray light curves of several blazars useful to study their variability on time scales from day to several 

37 months (see e.g. Abdo et al. 2010c). Daily light curves can be obtained for several blazars and for those 

38 exceptionally bright it is possible to observe significant occasional variations on timescales shorter than a 

39 day. We have now observed such events in 3C273, the nearest quasar. 

40 3C 273, the first quasar discovered by Schimdt (1963) and the first extragalactic source detected by 

41 COS-B in the gamma-ray band (Swanenburg et al. 1978), is one of the most extensively studied AGN 

42 across the entire electromagnetic spectrum. It is classified as a flat spectrum radio quasar (FSRQ) and has 

43 a redshift z = 0.158 (see e.g., Strauss et al., 1992). It was observed by EGRET (3EG J1229-I-0210 in 

44 Hartman et al. 1999) at an average flux of 0.18 x 10~^ ph cm~^ s~^with a peak flux of 1.27x 10~^ ph cm~^ 

45 s~^(Nandikotkur et al. 2007). Despite the huge amount of data collected to now (see for instance Soldi et 

46 al. 2008), its behaviour is still surprising and raises new challenges for physical models. It was detected by 

47 the LAT experiment since the beginning of the observation in 2008 (Marelli, 2008) and is identified with 

48 the gamma-ray source IFGL J1229.1-I-0203 in the First Fermi LAT catalog (Abdo et al. 2010a). For about 

49 one year its behavior was characterized by a rather stable flux of ~0.3 x 10~^ ph cm~^ s~^, with some flares 

50 superposed reaching peak level higher by about an order of magnitude. 

51 Beginning at the end of July 2009, 3C 273 started a brightening phase during which some very bright 

52 outbursts were observed. The first event lasted about 10 days in August 2009 (Bastieri, 2009) and was 
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53 characterized by fast rise and decay times (< 1 day) and a relatively stable 'plateau'. In contrast the second 

54 and third otbursts were sharply peaked from September 15 to 17 and from September 20 to 23, 2009 (Hill, 

55 2009), when 3C 273 reached a peak photon flux above 10^^ ph cm^^ s^^. Before the very recent (December 

56 2009) flare observed in 3C 454.3 (Escande & Tanaka, 2009; Stirani et al., 2009), it was, therefore, the 

57 brightest extragalactic source, non Gamma-Ray Burst, observed by Fermi. EGRET detected only a few 

58 blazars (e.g. 3C 279 and PKS 1622-297) with a flux above 10"^ ph cm"^ ^-i 

59 Thanks to the high flux is was possible to obtain light curves with a good S/N ratio with a time bin- 

60 ning of only six hours (corresponding to two scans of sky) and therefore we can describe the evolution of 

61 this outburst with a level of detail never reached before. In this letter we present the results of the LAT 

62 observations. 



63 2. Observations with the Large Area Telescope 

64 The Fermi-hAI is a pair-conversion 7-ray telescope sensitive to photon energies from 20 MeV to 

65 > 300 GeV. It is made of a tracker (composed of two sections, front and back, with different capabilities), 

66 a calorimeter, and an anticoincidence system to reject the charged-particle background. The LAT has a large 

67 peak effective area (~ 8000 cm^ for 1 GeV photons in the event class considered here), viewing ~2.4 sr 

68 of the full sky with excellent angular resolution (68% containment radius better than ~ 1° at E' = 1 GeV). 

69 A full description of the LAT instrument and its performance are reported in At wood et al. (2009). During 

70 the first year, the telescope operated in sky-survey mode, observing the whole sky every 3 hours. 

71 Only photons in "diffuse" class with energies greater than 100 MeV were considered in this analysis. 

72 We kept only events with reconstructed zenith angle < 105° in order to reduce the bright 7-ray albedo 

73 from the Earth. In addition, we excluded the time intervals when the rocking angle was more than 52° and 

74 when the Fermi satellite was within the South Atlantic Anomaly. The standard Fermi-LAT ScienceTools 

75 software packag^ (version v9rl5p6 ) was used with the P6 V3 set of instrument response functions. This 

76 version of the Science Tools takes into account the correction for the residual acceptance variation with 

77 count rate. Light curves were produced in 1-day (see Figure 1) bins over the period July-September 2009, 

78 and in 12 and 6-hour (see Figure 2) bins in September 2009. For each time bin, the flux, photon index 

79 and test statistic of 3C 273 were determined, in the energy range 0.1-100 GeV. We analyzed a Region of 

80 Interest (Rol) of 12° in radius, centered at the position of the 7-ray source associated with 3C 273 , using 

81 the maximum-likelihood algorithm implemented in gtlike. In the Rol analysis the sources were modeled as 

82 simple power-law {F = KE~^). The source model parameters are free for all point sources within 5° of 

83 3C 273 (extracted from the LAT IFGL catalog see Abdo et al. 2010a). In the Rol model, we also included 

84 all sources from 5° to 17° of 3C 273, with their model parameters fixed to their catalog values. 

85 The Galactic diffuse background model included in the Rol analysis is the currently recommended ver- 
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Fig. 1. — The 1-day gamma-ray light curve of 3C 273 in the period from July 1 (MJD=55013) to September 
30, 2009 (MJD=55104). 
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86 sion (glLiem_v02), publicly released through the Fermi Science Support Center. The isotropic background 

87 (including the 7-ray diffuse and residual instrumental backgrounds) model was derived from an overall fit 

88 of the diffuse component over the \b\ > 30° sky. A detailed documentation of the Galactic diffuse model 

89 and corresponding isotropic spectrum is available from the Fermi Science Support Centei[§. 

90 All errors reported in the figures or quoted in the text are 1-cj statistical errors. The estimated systematic 

91 uncertainty on the flux is 10% at 100 MeV, 5% at 500 MeV and 20% at 10 GeV. 



92 3. Time evolution of tlie outbursts 

93 The light curve of the very large outburst for photon energies above 0. 1 GeV and with a time binning 

94 of 6 hours is plotted in the upper panel of Fig. 2. The outburst's structure consists of two main flares, each 

95 one having a total duration of about a couple of days and peak flux of ~1.2x 10^^ ph cm^^ s^^. The first 

96 flare shows a nice smooth time profile, suggesting that it is well resolved; the second flare, after a first short 

97 peak, shows a series of high bins that can be indication either of substructures or of a long decay tail. 

98 The profile of the first flare is clearly asymmetric with a fast rise of the duration of about 12 hours, and 

99 a longer decay of about 36 hours. We modelled its evolution by means of the law of the following formula: 

F{t) = Fo + K{l + t/Tir{l-t/T2f (1) 

100 that is derived from the statistical Beta distribution (see, for instance Johnson, Kotz, & Balakrishnan 1994) 

101 and is also known as a Pearson curve of Type I (Smart 1958). It is defined in the interval (— Ti, T2) (a shift 

102 of the time origin is necessary for the fitting to data) and F{—Ti) = F(T2) = Fq at the extremes of that 

103 interval, then it is well suitable to model burst of finite durations. This function can be used to model the 

104 time evolution of flares, having either convex or concave rising and decaying profiles. It is more general 

105 than the one applied in Abdo et al. (2010c), which presented and analyzed the gamma-ray light curves of 

106 the 106 high-confidence Fermi LAT Bright AGN Sample (Abdo et al. 2009) obtained during the first 11 

107 months of the Fermi survey. The flares of the 10 brightest blazars were modelled by a law, defined as the 

108 inverse of the sum of two exponentials, to quantify the symmetry and duration of individual episodes. The 

109 total duration of the flare is given by Ti + T2, while we can easily measure its asymmetry by evaluating the 

110 difference between the decay and the rise time divided by their sum, which gives: 

a + p 

111 For the first burst we find Ti = 18 hours, T2 = 60 hours, a = 0.5, /3 = 3.2 and ^ = 0.73. The complex 

112 structure of the second flare make the fitting of Eq. (1) very hard, however, considering only the first five 

113 data points we obtained a reasonable fit with Ti = 12 hours and T2 = 33 hours and using the same values 



114 



of a and /? (see Fig. 2). The resulting x are smaller than unity in both cases. This result differs from those 
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115 of the other six and much longer bursts observed by the LAT to July 2009 (Abdo et al. 2010c). These flares 

116 had durations in the range from 8 to 58 days and were generally nearly symmetric, all they had negative 

117 values ^, usually in the interval (—0.38, —0.22) with the only exception being the longest one for which ^ 

118 was —0.72. Although the function used in the flare analysis of 3C 273 presented by Abdo et al. (2010c) was 

119 not that of Eq. (1), the quantity ^ was the same and depends only on the time elapsed between the peak and 

120 the starting and ending points of the flare. The highest brightness level reached in those flares, however, was 

121 much lower than that of the September 2009 outbursts: the typical flux of 3C 273 was around 3x10^"^ ph 

122 cm~^ s"^ with peak values during the flares of 2.5 x 10"^ ph cm~^ s~^. Comparable fluxes were measured 

123 by the GRID instrument, on board the AGILE satellite, from December 2007 to January 2008 (Pacciani et 

124 al. 2009). 

125 It is interesting to evaluate the rising and decay flux rates of the two flares. The rises were both very 

126 short, and the estimates of the corresponding flux rates at about half of the rising branch are 1.3x10-8 ph 

127 cm~^ s^^ h^^and 2.0 x 10^^ ph cm^^ s^^ h~^for the first and second flare, respectively. The decays evolved 

128 on longer time scale with flux rates of 0.3 xlO^^ ph cm^^ s^^ h^^and 0.5 xlO^^ ph cm^^ s^^ h^^. 

129 The very fast variability, however, does not correspond to an increase of the time-integrated energy 

130 output from the source, at least when they are compared with that of August. In fact the fluence of this flare, 

131 estimated by the integration of the best fit model over the duration of the event, is 2.06 ph cm~^, while those 

132 of the other and higher two flares are 1.59 and 1.05 ph cm^^. The last value is that of the main peak and 

133 increases by about 0.5 ph cm~^when the three last bins are included. 

134 A comparison of the power density spectrum (PDS) of the 90 day flaring episode in Fig. 1 with one 

135 computed for the first 11 months of data, did not show any substantial difference in shape, although the 

136 flare PDS had slightly more power at frequencies >0.1 days^^. The fractional variability (rms/flux) for 

137 frequencies above 0.01 days^^ in the PDS, was however, three times higher during the flaring episode than 

138 during the initial 11 months (0.99 and 0.32 respectivelly). In other words, the source was more variable 

139 during the flaring period of July-September than in the preceding 11 months. 

140 4. Spectral behavior 

141 Figured shows the spectral energy distribution, in the band 100 MeV - 20 GeV, of 3C 273 during the 

142 September outbursts. The spectra were obtained fitting a power law model over seven energy bins, with 

143 the spectral index kept fixed at the values obtained by the fit in the entire energy range over the same time 

144 interval. 

145 We also computed the spectra of 3C 273 integrating over the two weeks from September 1 to 13, during 

146 the quiescent phase preceeding the September events, and found that it is well described by a power law with 

147 photon index r=2.86±0.02 (see Figure O. The spectra of the two major outbursts integrated over the time 

148 intervals from 14 to 18 September and from 19 to 22 September 2009, showed flatter slopes r= 2.43±0.01 

149 and 2.52ib0.01, respectively, confirming that the average spectrum during the oubursts was significantly but 

150 moderately harder than in the quiescent state. As reported in Abdo et al.(2010d) the spectrum of 3C273 
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Fig. 2. — The upper panel shows the 6-hour light curve of 3C 273 during the large outbursts in September 
2009. The points marked with triangles indicate upper limits. The middle and lower panels show 12-hour 
light curves of 3C 273 during the large outburst in September 2009 in the energy range 0.1- 0.4 GeV and at 
energies higher than 0.4 GeV, respectively. 
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151 integrated over the 6-month showed a significant break around 1.5 GeV. Our spectra of the source integrated 

152 over shorter periods of time do not provide indication for a curvature in the spectra. 

153 It is to note that so far the spectral change with respect to flux has been observed to be very moderate 

154 in other FSRQs (e.g. 3C 454.3, see Abdo et al. 2010d), even during large flares. An interesting indication 

155 of a possible spectral evolution during the major bursts of 3C 273 can be obtained by the light curves in two 

156 energy bands. The central and lower panel of Figure 2 shows the outburst evolution in the energy ranges 

157 (0.1 - 0.4) GeV and above 0.4 GeV with a time bin width of 12-hour. Despite the large differences in 

158 the photon fluxes, we see that the profile of the peaks are rather similar with an indication that above 0.4 

159 GeV peaks are sharper and with shorter decay times. 

160 Time resolved spectral indices show possible spectral changes during the two major bursts as illustrated 

161 by the T vs Flux plots in Figure 4. The start and end ranges of the data used to demarcate the two flares were 

162 defined as the time intervals containing the 90% of the total flux, starting when the 10% level was measured. 

163 A "counterclockwise" loop is observed during the first flare. The spectrum becomes steeper in the declining 

164 phase and harder in the brightening phase. During the second flare the spectral index changes following a 

165 path that is essentially "clockwise" . In this case the brightening phase was very rapid and the spectral index 

166 remained constant, after the flare peak the spectrum becomes harder and then softer again. A clockwise loop 

167 can indicate that during the first flare the flux started to increase at low energy and then propagate to high 

168 energy, the inverse could have happened during the second flare. 

169 5. Discussion 

170 The large outburst of 3C 273 observed by Fermi in September 2009 revealed that blazars can reach 

171 very high brightness levels for quite short short time intervals of 1-10 days. 

172 The fluences of these events, however, are comparable to that of longer but less intense flares. We 

173 measured significant variations of the flux over time scales as short as six hours, which occurred with only 

174 a mild change of the spectral shape. 

175 In contrast to previous low intensity flares (Abdo et al., 2010c), the very strong events of September 

176 have light curves characterized by decay times longer than the rise time. Light curves at different energies 

177 provide some indication, in particular for the second flare, that the decay rate above 0.4 GeV is shorter than 

178 below. The most natural interpretation of this difference is that it is due to the radiative coohng of the high 

179 energy electrons responsible for the 7-ray emission. We can therefore to estimate that this time, as observed 

180 in the Earth's frame, is of the order of 0.5-1 day. Such short lifetimes explain why longer flares are generally 

181 symmetric. In fact, it is possible that they are structured in a series of short subflares with typical durations 

182 of two-three days, much shorter than the total flare length. The apparent rise and decay times of long events 

183 are then likely due to the superposition of quite shorter flares and any information on the radiative lifetimes 

184 is practically lost. 

185 In September 2009 the apparent position of 3C 273 was rather close to the Sun, making it was prac- 
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Fig. 3. — The Spectral Energy Distributions of 3C 273 in the quiescent phase and during the two September 
outbursts. 
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186 tically impossible to perform observations in other frequency ranges. Therefore we do not have data to 

187 Study the evolution of the broad band Spectral Energy Distribution. The typical SEDs of FSRQs (Abdo et 

188 al. 2010b) peak at rather low 7-ray energies, around 0.1 GeV and frequently lower. In the case of 3C 273, 

189 which was the most significant AGN detected by COMPTEL in the 1-30 MeV range, the peak of the inverse 

190 Compton component was in the range 1-10 MeV (Collmar et al 2000). Our data show that the photon index 

191 was always steeper than 2, and this indicates that the peak energy remained below 100 MeV also during the 

192 flares. 

193 7-ray activity has been found to be related to changes of the blazars' radio structure observed with 

194 VLBI. In this respect 3C 273 is one of the most interesting sources because of the flux and of the low 

195 redshift that allows a fine spatial resolution. Jorstad et al. (2001), on the basis of a large data set at 22 

196 and 43 GHz, found an association between the ejection of superluminal radio knots and high states of 7- 

197 ray luminosity in ten blazars in EGRET observations, including 3C 273. They concluded that both the 

198 radio and high energy events are originating from the same shocked region of a relativistic jet. Similar 

199 correspondences were already reported for the two much more distant FSRQs S5 0836+710 (Otterbein et al. 

200 1998) and PKS 0528+134 (Britzen et al. 1999). 

201 More recent results on the MOJAVE VLBA sample and Fermi-LAT observations (Lister et al. 2009, 

202 Savolainen et al. 2009a) provided evidence that 7-ray loud blazars have a Doppler factor higher than non 

203 LAT-detected sources. A detailed plot of the kinematics of the various components in the radio jet of 

204 3C 273 can be found in Lister et al. (2009b): the resulting mean superluminal velocity is ^app=13.4 with an 

205 estimated Doppler factor (5=16.8. 

206 It will be very interesting to verify if the exceptional outbursts of September 2009 will or will not be 

207 associated with the ejection of new superluminal knots, possibly with even higher velocity and Doppler 

208 factor. Moreover, the discovery of a possible connection between the peak intensity and rise time of the 

209 7-ray outbursts with the VLBI parameters can be very useful to constrain the modelling and the energetics 

210 of perturbations in the jet. 
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